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Abstract 

The quasigeostrophic equation is a prototypical geophysical fluid 
model. In this paper, we consider time-periodic motions of this model 
under dissipation and time-dependent wind forcing. We show that when 
the wind forcing is time-periodic and the spatial square-integral of the 
wind forcing is bounded in time, the full nonlinear quasigeostrophic 
model has time-periodic motions, under some conditions on f3 parame- 
ter, Ekman number, viscousity and the domain size. 
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1 Introduction 



The two dimensional barotropic quasigeostrophic (QG) equation has been de- 
rived as an approximation of the rotating shallow water equations by the 
conventional asymptotic expansion in small Rossby number ([11]). The low- 
est order approximation gives the barotropic QG equation, which is also the 
conservation law for the zero-th order potential vorticity. Warn et al. [It] and 
Vallis jL5[ emphasize that this asymptotic expansion is generally secular for 
all but the simplest flows and propose a modified asymptotic method, which 
involves expanding only the fast modes. The barotropic QG equation also 
emerges at the lowest order in this modified expansion. 

Schochet (Hl3"l) has recently shown that quasigeostrophy is a valid approx- 
imation of the rotating shallow water equations in the limit of zero Rossby 
number, i.e., at asymptotically high rotation rate. For related work in the 
three dimensional baroclinic QG model, see, for example, || , || , and @ . 

It is known that the linearized two dimensional barotropic QG equation 
without forcing and without dissipation has time-periodic solutions, and these 
time-periodic solutions are damped away by Ekman or viscous dissipation 
(0, pages 147 and 236). 

In this paper, we show that when the wind forcing is periodic and when 
its spatial square-integral is bounded in time, the full nonlinear QG equation 
with Ekman and viscous dissipation has time-periodic solutions. We use a 
topological technique from nonlinear global analysis (||10||). 



2 Quasigeostrophy 

The two dimensional barotropic QG equation is page 234) 

Ai/j t + J(V>, At/0 + (3ifj x = vA 2 if> - rAif; + f(x, y, t), (1) 

where ifj(x,y,t) is the stream function, (3 > is the meridional gradient of 
the Coriolis parameter, v > is the viscous dissipation constant and r > is 
the Ekman dissipation constant and f(x, y, t) is the wind forcing. Moreover, 
A = d xx + d yy is the Laplacian operator in the x, y plane and J(f, g) = 
fx9y — fy9x is the Jacobi operator. This equation is written in the often-studied 
case of infinite Rossby deformation radius and with flat bottom, applicable 
for planetary-scale solutions. The situation with infinite Rossby deformation 
radius is equivalent to the rigid-lid approximation (||, p223). 
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In this paper, we assume that the wind forcing f(x, y, t) is periodic in time 
with period T > 0. 

Introducing the relative vorticity u(x,y,t) = Aip(x,y,t), the equation (fl|) 
can be written as 

LO t + J (if), u) + /3i[) x = vAuj -ru + f(x, y, t) , (2) 

where (x, y) G D, an arbitrary bounded planar domain with piecewise smooth 
boundary. This equation is supplemented by zero Dirichlet boundary condi- 
tions for both tp and u = Ai/j, together with an appropriate initial condition, 
i.e., we require 

ip(x,y,t) = on dD , (3) 
u(x,y,t) = ondD, (4) 
u(x,y,0) = u (x,y) . (5) 

These boundary conditions have been used in analytical and numerical study 
of this model in, e.g., j|, W^M- 

We note that the Poincare inequality (M) 

[ g 2 (x,y)dxdy < - — - f \Vg\ 2 dxdy (6) 

JD 71 J D 

holds with these boundary conditions, where \D\ is the area of the domain D. 
The global well-posedness (smooth solutions) of this dissipative model can be 
obtained similarly as in, for example, []J, [|17| and [|1^ . 



It is well-known that the linearized two dimensional barotropic QG equa- 
tion without forcing and without dissipation 

A^ + /^ = 0, x, ye [0,1], (7) 

has time-periodic solutions p. 146-149). For example, assume the basin 

is the unit square on the (3— plane, with boundary condition ip — 0. The 
boundary is a streamline in this case. By separating variables, one can find 
that the equation ([7]) has time-periodic solutions (||11||, p. 146-149) 

ip mn (x,y,t) = cos(- h a mn t)sin(mTTx)sin(niTy),m,n = 1,2,3, ••• , (8) 

which are basin-scale traveling-wave oscillations with dispersion relation, 

7TITL I ' ^ ) 

2u\/m 2 + n 2 
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These are basin-scale normal modes (planetary waves) for QG in the rigid-lid 
approximation. Each mode consists of a carrier wave cos( ^ x +cr mn t) moving 
to the left (westward) and modulated by an envelope of sine functions which 
maintain the boundary conditions. These time-periodic linear solutions are 
damped away by Ekman or viscous dissipation ( Hill , p. 236). 

We address the question of whether there are any basin-scale time-periodic 
solutions in the full nonlinear dissipative QG dynamics (HD with time-periodic 
wind forcing. 



3 Dissipativity 

In the following we use the abbreviations L 2 = L 2 (D), Hq = Hq(D), . . . for the 
standard Sobolev spaces. Furthermore, let <•,•>, || • || denote the standard 
scalar product and norm in L 2 , respectively. We need the following inequalities 

(0)- 

Young inequality: 

AB<-A 2 + —B 2 , (10) 
~ 2 2e V ' 

where A, B are non-negative real numbers and e > 0. 

Gronwall inequality: If an integrable function y(t) satisfies that 

^<Ay + B, t>0, (11) 
for some constants A, B with A ^ 0, then 

y(t) < [?/(0) + j]e At - j, t>0. (12) 

We now show that the system (|2|), under boundary conditions (|3|), (f|), 
is a dissipative system in the sense (]TU| or 0) that all solutions u(x,y,t) 
approach a bounded set in L 2 (D) as time goes to infinity. A T— time-periodic 
dissipative system in a Banach space has at least one T— time-periodic solution. 
This result follows from a Leray-Schauder topological degree argument and the 
Browder's principle (0], p. 235). 

Multipling (0) by u> and integrating over D, we get 

_ ^ _ ll a, l| 2 + / J(ip,uj)udxdy + (3 ( ip x ujdxdy 
2 dt jd jd 

= —v \ I VcjI 2 — rllo'll 2 + / f(x,y,t)udxdy. (13) 

JD JD 
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Note that 

/ J(ip, u))udxdy = 0, (14) 
jd 

via integration by parts; see also |] or ||. Moreover, using the Young and 
Poincare inequalities, we get 



\(3 f ^ x udxdy\ < h{( #+/ u 2 ) 

JD 2 JD JD 

1 \D\ f f 

< ~/3( / uj 2 dxdy + / uj 2 dxdy) 

2 7T JD JD 



\pV— + 1) / w 2 dxdy, (15) 

2 7T JD 



— i/ / |Vu;| 2 c/x(iy < — - — - / u 2 dxdy. (16) 
JD \D\ JD 

We further assume that the square-integral of the wind forcing f(x,y,t) 
with respect to x, y is bounded in time. 

/ f(x,y,t)cudxdy<- [ f 2 (x,y,t)dxdy + e [ oo 2 dxdy < M + e / u 2 dxdy, (17) 
jd e Jd jd jd 

where e > is to be determined, and M > is a time-independent constant 
since we have assume that the square-integral of f(x, y, t) with respect to x, y 
is bounded in time. 

Putting O, (PD, (ELD, AHD into flT|), we obtain 



--T-IMI 2 + ol I io 2 dxdy<M, 
2 dt Jd 



where 



Assume that 



tip 1 ,\D\ 
a = i r+ m\-of 3 (— + 1 )-^ 

\JJ\ 2 TV 



\D\ 2 7T 

We can then take e > small enough such that a > 0. Thus, by the Gronwall 
inequality, we have 

lk|| 2 <(||o;o|| 2 )e- 2at + — . 20) 

a a 



Hence all solutions u> enter a bounded set {u : \\u)\\ < y 77} as time goes 
to infinity. The system (§) is therefore a dissipative system and hence has at 
least one T— time-periodic solution. 
We then have the following result. 

Theorem 1 Assume that 

(i) the wind forcing f(x,y,t) is time-periodic with period T > 0, and its 
square-integral with respect to x,y is bounded in time; and 

(ii) r + pjr > + 1), where (3 > is the meridional gradient of the 
Coriolis parameter, is > is the viscous dissipation constant, r > is the 
Ekman dissipation constant, and \D\ is the area of the bounded domain D. 

Then the dissipative quasigeostrophic model 

A^ t + J(iP,AiP)+/3^ x = uA 2 ^-rA^ + f(x,y,t) , (21) 

UJ = Ai/j , (22) 

ip(x,y,t) = ondD, (23) 

tv(x,y,t) = ondD, (24) 

u;(x,y,t + T) = u(x,y,t) , (25) 

has at least one time-periodic solution with period T > 0. 

We remark that it is generally difficult to show existence of time-periodic 
motions for a spatially extended evolution system. Our result provides such 
a proof of existence, under some conditions on (3 parameter, Ekman number, 
viscousity and the domain size, for a prototypical geophysical fluid model. 
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